arXiv:1508.00672vl [astro-ph.SR] 4 Aug 2015 


EVOLUTION AND DYNAMICS OF A SOLAR 


ACTIVE PROMINENCE 

T. Magara 

Dept, of Astronomy and Space Science, School of Space Research, Kyung Hee University, 
1732 Deogyeong-daero, Giheung-gu, Yongin, Gyeonggi-do, 446-701, Republic of Korea 

*E-mail: magara@khu.ac.kr 
Received ; Accepted 

Abstract 

The life of a solar active prominence, one of the most remarkable objects on the Sun, is 
full of dynamics; after first appearing on the Sun the prominence continuously evolves with 
various internal motions and eventually produces a global eruption toward the interplane¬ 
tary space. Here we report that the whole life of an active prominence is successfully re¬ 
produced by performing as long-term a magnetohydrodynamic simulation of a magnetized 
prominence plasma as was ever done. The simulation reveals underlying dynamic processes 
that give rise to observed properties of an active prominence: invisible subsurface flows self- 
consistently produce the cancellation of magnetic flux observed at the photosphere, while 
observed and somewhat counterintuitive strong upflows are driven against gravity by en¬ 
hanced gas pressure gradient force along a magnetic field line locally standing vertical. The 
most highlighted dynamic event, transition into an eruptive phase, occurs as a natural con¬ 
sequence of the self-consistent evolution of a prominence plasma interacting with a magnetic 
field, which is obtained by seamlessly reproducing dynamic processes involved in the forma¬ 
tion and eruption of an active prominence. A manuscript with high-resolution figures is found 
at http://web.khu.ac.kr/~magara/index.html. Here is a tilde. 
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1 Introduction 


Solar prominences are, just as their name tells, one of the prominent objects forming in the solar 
corona (Tandberg-Hanssen 1995). They are also called filaments when observed as a dark thin and 
long object on a solar disk. Prominences have various sizes with different appearances; their physical 
nature has been investigated in detail via theoretical modeling and multi-wavelength observations, 
which has revealed that a magnetic field frozen in a prominence plasma plays the role of a skeleton in 
determining the structural and evolutionary properties of a prominence (Kippenhahn & Schliiter 1957; 
Kuperus & Raadu 1974; Schmieder, Malherbe, Mein & Tandberg-Hanssen 1984; Martin, Livi & 
Wang 1985; van Ballegooijen & Martens 1989; Choe & Lee 1992; Low 1996; Aulanier & Demoulin 
1998; Martin 1998; DeVore & Antiochos 2000; Martens & Zwaan 2001; Chae, et al. 2001; Pevtsov, 
Balasubramaniam & Rogers 2003; Lites 2005; Mackay & van Ballegooijen 2005; Magara 2007; 
Okamoto et al. 2008; Berger et al. 2011; Priest 2014 and references therein). 

Prominences are classified into quiescent prominences and active prominences', the former are 
found outside solar active regions and size and lifetime are larger and longer than the latter which form 
in active regions and frequently produce explosive phenomena such as solar flares. This suggests that 
the formation processes of these two types of prominences may be different (Priest 2014). Here we 
report that dynamic processes giving rise to observed properties of an active prominence are seam¬ 
lessly reproduced by performing as long-term a magnetohydrodynamic simulation of a magnetized 
prominence plasma as was ever done. The simulation demonstrates how these dynamic processes are 
self-consistently driven and maintained. 

2 Simulation setup 

The simulation was performed in Cartesian coordinates ( x,y,z ) with the 2 -ax is directed upward. 
2 < 0, 2 = 0 and 2 > 0 correspond to the subphotosphere, the solar surface and the solar atmosphere 
extending above the surface, respectively. The simulation box is (-200,-200,-10) < ( x,y,z ) < 
(200,200,190) which is discretized into non-uniform grid cells whose size is (Ax, Ay, Az) = 
(0.2,0.2,0.2) for (-20,-20,-10) < ( x,y,z ) < (20,20,20) and it gradually increases up to (4,4,4) as 
|x|, \y\ and 2 increase. The total number of cells is N x x N y x N z = 383 x 383 x 242. The normaliza¬ 
tion units of the simulation are given by 2 H ph (length), c p h (velocity), p ph (gas density), p p hC 2 h (gas 
pressure), T ph (temperature) and (p p hC 2 h ) 1 / 2 (magnetic field), where H ph , c p h, p p h and T ph represent 
pressure scale height, sound speed, gas density and temperature at the photosphere. As the initial 
state of the simulation we adopt the same state as described in Magara (2012): a cylindrical magnetic 
flux tube of a Gold-Hoyle profile with a radius 77 = 2 and twist parameter b = 1 (Gold & Hoyle 1960) 
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is placed horizontally in the subphotosphere (the axis of the flux tube is set 4 below the y-axis). The 
flux tube keeps mechanical equilibrium with a background atmosphere stratified under the solar sur¬ 
face gravity g Q = 274 m/s 2 and a prescribed temperature profile based on a solar atmospheric model 
(Vernazza, Avrett & Loeser 1981). Boundary conditions and wave damping zones placed near all the 
boundaries are the same as the ones used in Magara (2012). The simulation is based on viscous mag¬ 
netohydrodynamic equations with the Reynolds number Re ~ 2.5 x 10 4 where a length scale is given 
by the typical height of an active prominence / ~ 7,000 km and a velocity scale is the free-fall veloc¬ 
ity v f = (2 g Q l) 1/2 20 km/s. The simulation started with the rise of the flux tube via the magnetic 

buoyancy (Parker 1955), which forms three emerging lobes of a magnetic field at intervals of 40, as 
shown by Figure la. Emergence of a twisted flux tube (Emonet & Moreno-Insertis 1998; Fan 2001; 
Magara & Longcope 2003; Archontis et al. 2004; Manchester IV et al. 2004; Galsgaard et al. 2005; 
Cheung & Isobe 2014 and references therein) is, observationally as well as theoretically, suggested 
to be one of the possible formation processes of an active prominence (Rust & Kumar 1994; Low 
1996; Magara 2007; Okamoto et al. 2008). The emerging portion of the flux tube that enters a low 
atmosphere (0 < z < 4) is subject to the Newton’s cooling characterized by the relaxation time r r ~ 1 
s (Stix 1991). 

3 Result and discussion 

Figures lb-d display snapshots of the middle emerging lobe taken at different evolutionary phases, 
which are viewed from the subphotosphere to see how a flow and magnetic field evolve below the 
surface when a prominence forms above the surface. These figures show a transition process from a 
divergent flow at an early phase during which the area of photospheric magnetic flux increases (figure 
lb) to a convergent flow at a late phase when positive and negative photospheric flux approach each 
other (figure Id), known as flux cancellation (Martin, Livi & Wang 1985). 

While flux cancellation is operating at the surface, a magnetic field emerging into the solar 
atmosphere develops a structure on the Sun, which is composed of the core field lines that form the 
main body of a prominence lying in the horizontal position and the envelop field lines overlying the 
main body (figure 2a). The gas density in the main body is about 10 5 or 10 6 times smaller than 
that in the photosphere, which is consistent with observations of active prominences (Priest 2014). 
Figures 2a-d show the evolution of selected field lines whose colors indicate the value of gas density. 
As the main body gradually changes shape from a horizontal shape into an fl shape, the gas density 
in the main body continues to decrease. This is because a prominence plasma easily falls down 
along an fl-shaped field line, which in turn enhances the magnetic buoyancy of fl-part of the field 
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line, facilitating the transition of field-line shape. The simulation demonstrates a positive feedback 
between the draining of mass and the transition of field-line shape, suggesting that it is inevitable for 
the prominence to proceed to an eruptive phase (figure 2d). 

Strong upflows of several tens of km/s exist in a prominence (Priest 2014) and their driving 
mechanism has been an issue to be clarified because the speed is comparable to the typical free-fall 
speed mentioned in section 2. In this simulation there are found locations of such strong upflows, one 
of which is indicated by the white box in figure 2d. Figure 3 presents a close-up view of this upflow 
location, where a field line is plotted to see the configuration of a magnetic field and distribution of 
gas pressure along the field line. The color gradient of the field line represents the variation of gas 
pressure, which explains that the upflows are driven by strong gas pressure gradient force along the 
field line locally standing vertical. 

The global eruption of an active prominence, which is an event at the final stage of evolution, 
is often associated with other explosive phenomena such as a solar flare and/or coronal mass ejection 
(Shibata & Magara 2011). They are transient but intense solar phenomena which potentially make a 
severe impact on the terrestrial environment. In order to understand a mechanism for producing these 
harmful phenomena and ultimately predict them, we have to clarify the dynamic properties of an 
active prominence evolving to global eruption. Seamlessly reproducing dynamic processes involved 
in the formation and eruption of an active prominence presented here gives an important step toward 
the comprehensive understanding and prediction of solar explosive phenomena. 

The author wishes to thank the Kyung Hee University for general support of this work. He 
also appreciates useful comments on the paper by S. Solanki. Visualization of simulation data was 
done using VAPOR (Clyne et al. 2007). This work was financially supported by the Basic Science 
Research Program (2013R1A1A2058705, PI: T. Magara) through the National Research Foundation 
of Korea (NRF) provided by the Ministry of Education, Science and Technology as well as the BK21 
program (20132015) through the NRF. 
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Fig. 1. (a): Snapshot of a local simulation domain (-50, -50,-10) < ( x,y,z) < (50,50,90) taken at t = 14 is presented. The green lines represent twisted 
field lines which form three emerging lobes while the gray-scale map at z = 0 gives the surface distribution of vertical magnetic flux in linear scale. Two color 
maps placed at x = 46 and y — 46 show the distributions of temperature and gas density in logarithmic scale. Gas density, temperature and magnetic field 
are normalized by photospheric density p p h, photospheric temperature T ph and (p P hc ph ) 1/2 where c ph is photospheric sound speed, (b)-(d): Snapshots 
of the middle emerging lobe taken at f = 20 (Figure 1b), 30 (Figure 1c) and 40 (Figure Id) show the evolution of selected field lines drawn in green. The 
gray-scale map placed at z = 0 shows the surface distribution of vertical magnetic flux in linear scale. The red arrows represent flow velocity field at y = 0. A 
unit vector scaled to c ph is given at top-left corner of each figure. 
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Fig. 2. (a): Snapshot of a local simulation domain (-80, -80, -10) < ( x,y, z) < (80,80,90) taken at t = 60.2 is presented. The colored lines represent 
selected field lines; colors indicate the value of gas density in logarithmic scale. The horizontal gray-scale map at z = 0 shows the surface distribution of 
vertical magnetic flux in linear scale. The vertical gray-scale map at x = 0 shows the distribution of gas density in logarithmic scale. The yellow arrows 
represent flow velocity field at z = o and a unit vector scaled to c ph is given at bottom-right corner. The blue arrows represent flow velocity field at x = 0 and 
a unit vector scaled to 10 c ph is given at top-right corner, (b)-(d): Snapshots of the same selected field lines as those in Figure 2a are presented, except for 
t = 70.2 (Figure 2b), f = 80.2 (Figure 2c) and t — 90 (Figure 2d). The white box in Figure 2d shows a location of strong upflows, the close-up view of which 
is presented in Figure 3. 
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Fig. 3. Close-up view of the upflow location indicated by the white box in Figure 2d is presented. The color gradient of the displayed field line shows the 
variation of gas pressure in linear scale, normalized by photospheric gas pressure P p h = Pphc^bj -1 where 7 is the specific heat ratio. A unit vector scaled 
to 10 c P h is given at bottom-left comer. 
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